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heavy fermions 

PACS 68.37. Ef - Scanning tunneling microscopy (including chemistry induced with STM) 

PACS 73.20.-r - Electron states at surfaces and interfaces 

Abstract - We study theoretically a square lattice of the organometallic Kondo adsorbate iron(II) 
phtalocyanine (FePc) deposited on top of Au(lll), motivated by recent scanning tunneling mi¬ 
croscopy experiments. We describe the system by means of an effective Hubbard-Anderson model, 
where each molecule has degenerate effective d— orbitals with xz and yz symmetry, which we solve 
for arbitrary occupation and arbitrary on-site repulsion U. To that end, we introduce a gener¬ 
alized slave-boson mean-field approximation (SBMFA) which correctly describes both the non¬ 
interacting limit (NIL) [7 = 0 and the strongly-interacting limit U —» 00 , where our formalism 
reproduces the correct value of the Kondo temperature for an isolated FePc molecule. Our results 
indicate that while the isolated molecule can be described by an SU(4) Anderson model in the 
Kondo regime, the case of the square lattice corresponds to the intermediate-valence regime, with 
a total occupation of nearly 1.65 holes in the FePc molecular orbitals. Our results have important 
implications for the physical interpretation of the experiment. 


The advances in nanotechnology and the degree of con¬ 
trol of different parameters in systems with strongly inter¬ 
acting electrons provide a rich and challenging field for the 
theoretical understanding of strongly-correlated materials, 
and for future applications in spintronics and quantum in¬ 
formation [l|. The strong repulsion between 3 d electrons 
in the 3d series of transition-metal atoms (TMAs) has led 
to the observation of the Kondo effect in several systems 
in which these atoms GHU or molecules containing them 
[SHHi are absorbed on noble-metal surfaces. In its simpler 
form, the Kondo effect arises when a TMA has a (nearly) 
integer occupation of an odd number of electrons and spin 
1/2, and this spin is “screened” by a cloud of conduc¬ 
tion electrons in the metal, building a many-body singlet 
ground state. As a consequence, a narrow Fano-Kondo 
antiresonance (FKA) appears in the differential conduc¬ 
tance in scanning tunneling spectroscopy (STS) experi¬ 


ments. More exotic Kondo effects involving even occupa¬ 
tion and spin 1 jHJG2|, and orbital degeneracy [HUITIITSI 
were also observed. In particular, it has been shown that 
STS for a system consisting on a single iron(II) phtalocya¬ 
nine (FePc) molecule (see the inset of Fig. 1) deposited 
on top of clean Au(lll) (in the most usual “on-top” con¬ 
figuration) can be explained in terms of an effective SU(4) 
impurity Anderson model m- The four degrees of free¬ 
dom correspond to the intrinsic spin 1/2 degeneracy times 
two-fold degeneracy between molecular orbitals with the 
same symmetry as the 3 d xz and 3 d yz states of Fe, where z 
is the direction perpendicular to the surface. The Kondo 
model is the odd-integer valence limit of the impurity An¬ 
derson model, in which localized electrons (like those of 
the 3d shell of a TMA) are hybridized with a conduction 
band. 

Experiments in which an array of effective “impurities” 
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(TMAs or molecules containing TMAs) are deposited on 
metal surfaces constitute a further challenge to the theory 
because they should be described by the periodic Ander¬ 
son model, for which neither exact solutions nor accurate 
treatments, like the numerical renormalization group, ex¬ 
ist. Tsukahara et al. studied the STS of a square lattice 
of FePc molecules on Au(lll) pjj. They observed a split¬ 
ting of the FKA near the Fermi level. A first hint to 
the correct description of the system came soon when the 
corresponding experiments for an isolated FePc molecule 
could be well described by the effective SU(4) impurity 
Anderson model m- a natural step to construct a lat¬ 
tice model is to introduce an effective hopping between 
the molecules, which leads to a two-band Hubbard model 
when the hybridization H mix with the metallic states is ne¬ 
glected. When H mix is included, one obtains an orbitally- 
degenerate Hubbard-Anderson model. Lobos et al. |T6] 
derived this model and solved it using a slave-boson mean- 
field approximation (SBMFA) for infinite on-site repulsion 
U [17] . These authors were able to reproduce the observed 
splitting of the FKA, which in simple terms can be inter¬ 
preted as arising from split van Hove singularities due to 
the anisotropic hopping terms in the Hubbard part of the 
model. 

However, there are several indications that the theory 
needs to be improved in order to reproduce the experi¬ 
mental details. For instance, the energy necessary to add 
the first hole Ed « —0.1 eV [see Eq. (1)] was estimated 
from ab-initio calculations in the local density approxi¬ 
mation (LDA), which obtained spectral density of Fe 3 d xz 
and 3 d yz states 0.1 eV above the Fermi energy ep, which 
we set as the origin of energies (ep = 0) However, as de¬ 
scribed below, more accurate LDA+U calculations esti¬ 
mate Ed ~ —1.6 eV. Since the value U = 1.8 eV was 
estimated [10], Ed + U is small and in principle double 
occupancy of holes cannot be neglected contradicting the 
assumptions in Ref. [TB] (the hole occupancy is between 
1 and 2). In addition, that theory does not allow to ob¬ 
tain the position of the split dips in the FKA at the same 
position as in the experiment, and an ad hoc voltage (hor¬ 
izontal) shift in the differential conductance dl/dV was 
necessary to fit the experiment. While a shift in dl/dV 
(vertical) might be expected due to uncertainties in the 
background contribution due to conduction electrons, a 
voltage shift is more difficult to justify. As a comparison, 
for a Co atom on Au(lll) a SBMFA of the SU(2) impurity 
Anderson model provides an excellent fit of the observed 
dl/dV without the need of a voltage shift [15 

The above discussion shows the need of accounting for 
a finite U and multiple occupancy at the impurity sites. 
Kotliar and Ruckenstein originally implemented this idea 
within the SBMFA m, which was later generalized for 
orbitally-degenerate systems [201 - I22] . and constitutes a 
very successful tool for describing highly-correlated sys¬ 
tems [251125] . However, as we show below, this approach 
does not lead to the correct Kondo temperature for an 
isolated molecule in the Kondo regime. 


In this work we study the orbitally degenerate Hubbard- 
Anderson model as an effective model to describe a square 
lattice of FePc molecules on Au(lll), for realistic param¬ 
eters derived from the LDA+U. We construct a SBMFA 
that gives the correct Kondo temperature for one FePc 
molecule, and also reproduces very accurately the non¬ 
interacting limits (NILs) U = 0, for both the impurity 
and the lattice cases. Although physically we expect U to 
be rather large (i.e., approx. 1.8 eV), reproducing the ex¬ 
actly solvable NIL constitutes an important “sanity check” 
for the entire formalism and is important when multi¬ 
ple occupancies cannot be neglected. To the best of our 
knowledge, this is a novel theoretical development which 
improves over previous SBMFAs for orbitally-degenerate 
correlated systems [2H1122] . We also stress that the formal¬ 
ism presented here is in principle generic and applicable 
beyond the case of FePc/Au(lll) molecules. 

Our method allows for quantitatively accurate fits of 
the experimental dl/dV , both for the isolated molecule 
and for the square lattice. The fitting procedure allows to 
extract the parameters of our model. From here we con¬ 
clude that while the single FePc molecule is in the Kondo 
limit, the lattice of FePc molecules is in the intermediate- 
valence regime. 

The Hamiltonian was derived in Ref. [16]. We refer the 
reader to that work for details. It can be split into three 
parts. H mo i is a degenerate Hubbard model that describes 
the effective molecular states and the hopping between 
them, H c contains the conduction states, and H m i X is the 
coupling between molecular and conduction states: 
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The operators h" a annihilate a hole (create an electron) 
in the state \v Tijta ), where v = (a :,y) denotes one of the 
two degenerate molecular states with spin a at site with 
position Tij = iai + j a 2 of the square lattice, with a^ 
the Bravais lattice vectors in the directions x and y re¬ 
spectively. For simplicity we denote as x and y the or¬ 
bitals with symmetry xz and yz. The operator for the 
total number of holes at the molecule lying at site i+ is 
+« =_Ecry n ry,o-> with ,* ■ Th e mean¬ 

ing of P is x = y, and y = x. The hopping £2 between x (y) 
orbitals in the x ( y ) direction is larger than the hopping 
ti between x ( y ) orbitals in the y (x) direction. 
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H c corresponds to a band of bulk and surface conduc¬ 
tion electrons of the substrate. The operator c Tij ^ tCr anni¬ 
hilates a conduction hole with spin a and quantum number 
£ at position r tJ . Note that the form of H mix and H c as¬ 
sumes that the molecular states of each Hubbard site r {j 
is hybridized with a different conduction band. This is an 
approximation which is valid if the distance between the 
Hubbard sites is R 1 /ftp, with kp the Fermi momen¬ 
tum of the metallic substrate [261428] . as it is the case in 
our problem. In other words, the molecular states are well 
separated and if the hopping between them is neglected, 
the system behaves as dilute system of impurities. 

We use the subscript p to denote any of the four states 
at each site (x |, x 4, y t, y I)- The basic idea of the SB 
approach is to enlarge the Fock space to include bosonic 
states which correspond to each state in the fermionic de¬ 
scription. For example, the vacuum state at site i is now 
represented as ej|0), where e\ is a bosonic operator corre¬ 
sponding to the empty site; similarly s| // |0) represents 
the simply occupied state with one hole with “color” 77, 
d\ V 2 fim I corresponds to a state with double oc¬ 
cupancy and similarly for triple and quadruple hole oc¬ 
cupancy (we denote the bosons as t\ and q respec¬ 
tively). The fermion operators entering Eq. (|T|! in the 
new representation are given by 




nj,v 
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and where 

we 

assume z Tijtv = z Xi . 

,7) for 


the moment . The bosonic operator zj. t] corresponds to 
the creation of the fermion p at site r,; 7 -, in the bosonic 
sector of the Hilbert space. 

To restrict the bosonic Hilbert space to the physical 
subspace, the following constraints must be satisfied at 
each lattice site (in what follows we drop the site index 
rij for simplicity) 


1 — e^e + s^s v + dj /iV2 d Vir/2 

7) 1)1 <7)2 
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The interaction term in Eq. © can be written as 
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The advantage of the SB representation is that using Eq. 
(fo|l in the saddle-point approximation, in which all bosonic 
variables are replaced by numbers, the problem is reduced 
to a non-interacting one, in which the values of the bosons 
are obtained minimizing the free energy. A shortcoming of 
this procedure is that using Eqs. (BED , the exactly solvable 
NIL U = 0 is not recovered after performing the mean-field 
approximation (MFA). In this limit, all bosonic numbers 
become independent of position r^ and color 77 , they can 
be chosen real and positive and for a total fermion occu¬ 
pation n = {n Ti .), calling p = n/4, one has e 2 = (1 — p) 4 
(probability of finding an empty site), s 2 = p(l — p) 3 , 
d 2 = p 2 { 1 — p) 2 , t 2 = p 3 ( 1 — p) and q 2 = p 4 . This leads 
to z 3 v = [p( 1 — p)] 1 ' 2 , while the correct result is z 3 = 1 
for U = 0 . To remedy this problem for the SU(2) case, 
Kotliar and Ruckenstein have replaced the z 3 operator 
for another one Zi v which is equivalent to z 3 in the re¬ 
stricted Hilbert space, but when evaluated in the MFA 
for U = 0 gives z v = 1 m They have shown that with 
this correction, the SBMFA is equivalent to the Gutzwiller 
approximation. This procedure was generalized for an ar¬ 
bitrary number M of colors [20U221 . For M = 4 we can 
write (dropping again the site indices for simplicity) 


z v = 
O v L = 


— a 2 

Or)R = 

— a 3 


n - 1/2 z o 0 -V 2 
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where in Eq. © all pi ^ p. We have written the coeffi¬ 
cients Aj in Eqs. (17151) so that the formalism is electron- 
hole symmetric, giving the same z v for occupation n and 
4 — n. Note that z 3 gives a non vanishing result on states 
which do not contain p, while all terms of Eq. © except 
the first do contain this particle. Therefore z 3 0~^ 2 = z 3 . 
Similarly, since z 3 creates this particle in the bosonic lan¬ 
guage, and all terms except the first in Eq. ([ 8 |) project 
over states which do not contain it, 0~^ 2 z 3 = z 3 . Then 
Zri = z 3 in an exact treatment. However this is no more 
true in the MFA. Choosing Aj = 1, the second of Eqs. 
dU) gives {O v ii) = 1 — p, and by electron-hole symmetry 
(O v l) — P- This implies (z v ) = 1 for U = 0 recovering the 
correct NIL in the SBMFA for any occupation. 

Nevertheless, the widely used choice Aj = 1 leads to an 
incorrect Kondo temperature Tk for the impurity case in 
the limit U — > + 00 . From exact Bethe ansatz results, the 
Kondo temperature in this limit for the SU(N) Anderson 
model is (except for a factor of the order of 1 that depends 
on the way in which it is defined) |29j 
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T*~£>exp(^0, (9) 

where a constant density of conduction states p per color 
(spin and symmetry) extending from —D to D and the 
Kondo regime —Ed A are assumed, with A = npV 2 . 
The problem is reduced to an effective one-fermion model, 
in which the hybridization term is reduced by the factor 
(zi V ) = z , the term Hjj is added, and the constraints of 
Eqs. QJ should be added. Using the first constraint one 
can eliminate one of the bosons, and the second one is 
taken into account by introducing a Lagrange multiplier 
A, leading to an effective impurity level at Ed = Ed + A. 
The Green function of the pseudofermions is 

Gf(w) = (co-Ed + iA)- 1 , (10) 

where A — z 2 A is one of the definitions used for Tk ■ The 
change in energy due to the impurity becomes [29] 
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— [—4A + 2A In 
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+ 4 Ed arctan 



—4A(s 2 + 3d 2 + 3 1 2 + q 2 ) + 6 U(d 2 + 21 2 + q 2 ).{ 11) 


For U —> + 00 , one has d = t = q = 0, and one can 
use e 2 = 1 — 4 s 2 [from Eq. (©)]. The condition dE/d A = 
0 gives Ed = Acot( 7 rs 2 ) and minimizing Eq. (fill) with 
respect to s one obtains 


2 Adz 2 
7 r ds 
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In the Kondo limit, —Ed A, the total occupation n —> 
1, and then s —> 1/2 and Eq. (fl2ll gives 


A = 


D nEd 

7i exp ^A’ 


dz 

= "ds" s=1 / 2 " 


(13) 


In order for this equation to be consistent with Eq. ©, 
minus the derivative of z 2 evaluated at s = 1/2 should be 
C = N = 4. Using Eqs. (1316171 and [EJ, and taking A 4 = 1 
(to have a finite z for n —> 0 ) one obtains 


compensates the factor yjp of zq, while O v l —> 1. Then, 
the condition (7 = 4 implies A\ = 0. To determine the 
remaining coefficient we impose that the NIL for n — 1 
is also recovered. This leads to a quadratic equation for 
A 2 . Its solution nearest to 1 gives A 2 = 1.3185. Note that 
because of electron-hole symmetry, z — 1 also for n = 3 
and n —> 4. We have represented the resulting z as a func¬ 
tion of n for U = 0. It is a very flat function with a value 
near 1 , and the maximum deviation is at n — 2 for which 
z = 0.989. Therefore, our SBMFA reproduces accurately 
not only the Kondo limit for one impurity at large (7, but 
also the NIL in the general case. Another possibility to 
optimize the Aj is to relax the correct NIL for n —> 0 
and n —> 4 and to ask that the formalism captures the 
correct Kondo temperature Tk — D exp[—7 t( 7/(16A)] for 
Ed = —3(7/2 such that the system is in the electron-hole 
symmetric case (with n = 2) for large but finite U. This 
condition would imply A 2 = 2 ± y/3/2. However, for re¬ 
alistic parameters (as described below) this estimate gives 
Tk - 5x 10“' K which is too small, indicating that these 
values of Ed are not realistic. Thus, we keep the values of 
Aj described above. 

The differential conductance G(I4) = dl/dVf, measured 
in STS experiments, is a non-equilibrium process in which 
the conduction states of the tip of the scanning-tunneling 
microscope at a bias voltage 14 hybridizes with a linear 
combination l Tij of molecular and conduction states. This 
amounts to adding a perturbation V t ry + H.c.) to 

the Hamiltonian, where tk refers to the tip states and 


= 


E„(Cr 


where a is a normalization 


'-'Tijt'q T qdriji’r)') 

constant. In the limit in which V) is small, G(I4) is propor¬ 
tional to the density of states pi(eVb) of the states l rij at 
the position r ,; 7 of the tip [18]I30] . Here q controls the ratio 
of the hybridization of the tip with the molecular states 
with respect to the corresponding value for the conduction 
states (this interference is the physical origin of the Fano 
effect). Once the Green function of the molecular states 
G£ ijV (u) = z r [see Eq. (HI) for the impurity 
case] is obtained from the SBMFA, pi{w) = —Im[G;(w)]/ 7 T 
is given by a simple expression obtained using equations 
of motion [TB] . One has 


G 


1 — 4s 2 

(4 — 3A 3 )(1 — A\s 2 ) ’ 
16 

(4 — 3A 3 )(4 — Ai)' 


U —y -f*oo, 


(14) 

(15) 


Thus, using Aj = 1 leads to an incorrect exponential de¬ 
pendence of Tk on Ed for U —> +oo. To remedy this, we 
modify the coefficients Aj so that the condition G = 4 
is satisfied, and at the same time we try to recover the 
NIL U = 0 as accurately as possible. To reproduce the 
NIL for n = 4p —> 0, in which e 2 ~ 1 — 4 p, s 2 ~ p , and 
the rest of the bosons can be neglected, we must choose 
A 4 = A 3 = 1. This gives O v r = p [see Eq. ©] and 


G l rij (w) — y|[G£,» + (VG^. v (uj) + q) 2 G{... iI 7 (w)], 

(16) 

where G° c _ is the Green function of the conduction states 
before adding the molecules (we take that corresponding 
to a constant density of states p = 1/(2 D) extending from 
— D to D). 

According to the estimations of LDA+U (see supple¬ 
mental material of Ref. [10]), for an isolated FePc molecule 
on Au(lll), U = 1.84 eV and the energy of the localized 
molecular orbital is Eq = —3.91 eV. Therefore, the energy 
necessary to add the first hole on a molecular orbital is 
Ed = — (Eq + 3(7) = —1.61 eV. In the rest of this work we 
take U = 1.8 eV and for this system we take Ed = —1.6 
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eV (for the square lattice of FePc molecules we modify Ed 
as described below). We have chosen D = 3.65 eV and 
p = 1/(2 D). However, the relevant parameter which con¬ 
trols the hybridization is the product A = irpV 2 , while the 
spectral density of the molecular states, as well as pi (w)/ p 0 
depend very weakly on D. 
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Fig. 1: change of the spectral density of the operator sensed by 
the STM tip l (see text) as a function of frequency for an iso¬ 
lated FePc molecule. The circles correspond to 3.1G(a>/eV) — 
2.5 where G in arbitrary units was taken from Fig. 3(a) of 
Ref. m- Parameters are U = 1.8 eV, Ed = —1.6 eV, A = 11 
meV and q = —0.03. The inset shows a scheme of the FePc 
molecule. 

To describe the case of a single FePc molecule, we take 
A and q as fitting parameters. The resulting pi{ui), except 
for an additive constant, is proportional to the observed 
G. In Fig. [Q we represent the change A pi(oj) in pi(co) af¬ 
ter addition of the impurity [last term in Eq. (fTBl) ] which 
should be proportional to the change in G, and is less sen¬ 
sitive to the details of the conduction band. From the fit 
we obtain A = 11 meV and q = —0.03. Note that the 
agreement with experiment near the Fermi energy is very 
good. It is certainly better than the comparison between 
two different experimental realizations [Figs. 3 (a) and (b) 
of Ref. TD] for zero magnetic field]. For the parameters 
of our fit, we obtain a total occupation of 1.01 holes (2.99 
electrons) in the molecular states. The total probability 
of double hole occupancy is 1.3% and that of having no 
holes (4 electrons) is 0.14%. Therefore, we conclude that 
the single FePc molecule is in the Kondo regime. These 
results agree with those of Minamitani et al. na who used 
the numerical-renormalization group to interpret the ex¬ 
periments. However, those authors did not attempt to fit 


the experiment. 

To model the square lattice of FePc molecules on 
Au(lll), we need to introduce the intermolecular hopping 
elements fj. We assume £ 2/^1 = 3 as estimated previously 
m- The SBMFA again reduces the problem to an effec¬ 
tive non-interacting one. The latter has the same form as 
that explained in the supplemental material of Ref. [16] . 
but in our case a finite U is used and therefore the quasi¬ 
particle weight z 2 .. rj is different and given by Eqs. mm 
and|U) with A\ = 0, A 2 = 1.3185, and A 3 = A4 = 1 as de¬ 
scribed above. Although for the lattice, the model is not 
SU(4) symmetric, the effective non-interacting problem re¬ 
tains this symmetry m and in the homogeneous SU(4) 
symmetric solution z Tij , v = z independent of site, spin and 
orbital. In order to compare with the experimental results 
for the square lattice [5], we retain the same values of U 
and A as for the impurity case, and use t\, q and Ed as fit¬ 
ting parameters, q depends on the distance from the STM 
tip to molecules and the surface and therefore one expects 
that it is different from the one-molecule case. Concerning 
Ed since the molecules are negatively charged as they are 
placed on the surface uniEa, one expects that Ed low¬ 
ers for the lattice due to interatomic repulsion, since holes 
are stabilized. This agrees qualitatively with our findings. 
The quantitative aspects will be discussed below. 

The resulting change in the spectral density sensed by 
the STM tip is shown in Fig. [2] From the fit we ob¬ 
tain t\ = 7.5 meV, q = —0.006 and Ed = —1.828 eV. 
The latter is a surprising result since Ed+ U = — 28 meV 
has become slightly negative, favoring double occupancy 
but the system is in the intermediate valence regime, be¬ 
cause this energy is of the order of A = 11 meV. This 
value of Ed is imposed by the position in voltage of the 
observed differential conductance. Higher values of Ed in¬ 
crease the electron occupation and shift the structure to 
the left. For the parameters of Fig. [2] we obtain a total 
occupancy of 1.65 holes (2.35 electrons). The total single 
hole occupancy is 0.34, the total probability of double hole 
occupancy is 0.65, and other states can be neglected. 

The drastic change in the charge of the molecular xz 
and yz orbitals is surprising. We should note that first- 
principle studies of the charges for isolated FePc urn 
and CoPc [55] on Au(lll) indicate that in addition to 
these orbitals, also the 3d orbitals of the TMAs with sym¬ 
metry 3 z 2 — r 2 are partially occupied. Nevertheless, Mi¬ 
namitani et al. showed that the spin of these orbitals is 
screened in a first-stage Kondo effect at higher energy 
and that an effective model containing only xz and yz 
orbitals as the localized ones describes the low-energy 
physics m- When an isolated FePc molecule is placed 
on the metal surface, a partial transfer of electrons takes 
place from the metallic states and also from the xz and 
yz orbitals to the 3 z 2 — r 2 [32]. Our results suggest that 
the latter charge transfer is enhanced for a lattice of FePc 
molecules on Au(lll). The intra-atomic inter-orbital re¬ 
pulsion should favor this procedure, since a larger charge 
in the 3 z 2 — r 2 orbital increases the energy for occupying 
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Fig. 2: same as Fig. 1 for a square lattice of FePc molecules. 
The circles correspond to 0.28 G(cj/eV) — 0.7 where G in arbi¬ 
trary units was taken from Ref. [8]. Parameters are t/ =1.8 eV, 
Ed = —1.828 eV, A = 11 meV ti = 7.5 meV, and q = —0.006. 

xz and yz orbitals. These arguments suggest that a more 
general model that includes 3 z 2 — r 2 localized states might 
explain the change in the occupancy of xz and yz orbitals. 

In summary, we have studied an array of FePc molecules 
deposited on top of Au(lll), described by an effective 
Hubbard-Anderson model with degenerate effective or¬ 
bitals with xz and yz symmetry. To that end, we have 
introduced a generalized slave-boson mean-field approxi¬ 
mation (SBMFA), which correctly describes both the non¬ 
interacting (U = 0) and the strongly-interacting (U —> oo) 
limits. This is an important improvement over previ¬ 
ous formulations of the SBMFA for strongly-correlated 
orbitally-degenerate systems, which fail to describe the 
Kondo limit. We stress that our method is generic and 
has applications beyond the present case of FePc/Au(lll). 
Our results indicate that while the isolated FePc molecule 
can be described by an SU(4) Anderson model in the 
Kondo regime, the case of the square lattice corresponds 
to the intermediate-valence regime, with a total occupa¬ 
tion of nearly 1.65 holes in the FePc molecular orbitals. 
This conclusion is imposed by the position in voltage of 
the double-dip structure observed in dl/dV and is inde¬ 
pendent of the details of the parameters. We believe that 
the shift in Ed is partly due to intermolecular repulsion, 
but it might also indicate a redistribution of charge among 
the Fe 3 d electrons that is beyond our effective model. 
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